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Failure Analyses of Composite Bolted Joints 


ABSTRACT 


The complex failure behavior exhibited by composite 
bolted joints of graphite-epoxy (Hercules AS/3501) has 
been investigated for the net tension, bearing and shearout 
failure modes using combined analytical and experimental 
techniques. Plane stress, linear-elastic, finite-element 
methods were employed to determine the two-dimensional 
state of stress resulting from a loaded hole in a finite- 
width, semi-infinite strip. The stresses predicted by 
the finite-element method were verified by experiment to 
lend credence to the analysis. Both numerical and 
experimental parametric studies were conducted to 
investigate the influence of joint geometric parameters 
on the state of stress and resultant strength of the 
joint. The resulting functional relationships found to 
exist between bolted joint strength and the geometric 
parameters, bolt diameter (D), e/D and W/D, were applied in 
the formulation of semi-empirical strength models for the 
basic failure modes. A point stress failure criterion 
was successfully applied as the failure criterion for 
the net tension and shearout failure modes. The bearing 
failure rode could not be accurately modeled using 
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two-dimensional plane-stress, finite-element analysis 
in combination with any known failure criterion although a 
maximum bearing stress criterion serves as a good order of 
magnitude estimate. Thus, using the results of the 
parametric studies, experimental test data and the point 
stress failure criterion, shearout and net tension 
strength models have been developed which account for the 
effects of geometry in composite bolted joint strengths. 
The implications of the results from this study on the 
further development of comprehensive bolted joint failure 
models for use in the efficient design composite bolted 
joints is discussed. 
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NOMENCLATURE 


Term 


Description 


edge distance (e) 

hole diameter (D) 
width (W) 

pin loading 

out-of-plane constraint 

bearing failure 
shear-out failure 


net tension failure 


ultrasonic "C" scan 


point stress failure 
criterion 


Distance from center of hole to 
the end edge of the laminate. 

The diameter of the fastener hole. 

Width of the coupon which is 
equivalent to the half spacing 
between fasteners. 

Reacting the joint loads with a 
pin which offers no out-of-piane 
restraint. 

Reacting the joint loads through 
the pin while constraining the 
out-of-plane deformations (defor- 
mations normal to plane of lami- 
nate) . 

A failure where localized crush- 
ing of the material reacting load 
leads to an elongation of the 
hole . 

Failure which s characterized by 
extensive cleavage damage parallel 
to the loading extending from the 
two sides of the fastener hole to 
the edge of the coupon. 


Normal tensile failure occurring 
in the material adjacent to the 
fastener hole. 

A nondestructive inspection tech- 
nique which uses ultrasonic waves 
to reveal information about the 
structure of the material. 

A failure criterion which predicts 
failure when the stress adjacent 
to a notch at a critical distance 
d 0 reaches the unnotched laminate 
strength. 
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INTRODUCTION 


Advanced composites are being considered for many 
load bearing structural applications requiring mechanically 
fastened joints. The cutouts required for mechanical 
fasteners always result in reduced structural load carrying 
efficiency and in composites this problem is amplified by 
the anisotropic nature of the material. The material 
anisotropy and brittleness intensify the stress concentration 
effects and provides for competing failure modes, each 
of which may be catastrophic. As a result the failure 
strength and mode of composite bolted joints are strongly 
dependent upon the material properties, laminate configura- 
tion and joint geometry. The development of a comprehensive 
composite bolted joint analysis procedure hinges upon 
understanding the influence of these parameters upon 
bolted joint strength and failure mode. 

The development of a comprehensive composite bolted 
joint analysis procedure involves the analytic determina- 
tion of the state of stress/strain in the joint element and 
application of an appropriate failure criterion. The closed- form 
solution to the elasticity problem of a loaded hole with 
a rigid inclusion in a finite-width, semi-infinite, 
anisotropic plate is intractable, hence numerical methods 
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must be used to determine the state of stress in the 
joint. Both finite-element [1-41 and approximate elasticity 
solutions (5,6) have been developed for the plane stress 
analysis of the stresses in an anisotropic plate loaded 
by a frictionless pin. The choice of an appropriate 
failure criterion has been the deficiency of most analyses. 
Distortional energy [1,2], Tsai-Wu i 3 ] , Maximum stress, 
Fracture Toughness [7] , and Average Stress failure criteria 
[8] have all been employed with various degrees of success. 
An acceptable degree of accuracy and reliability has not 
been obtained with any one of the above criteria for all 
failure modes primarily because there are different failure 
mechanisms associated with each mode and no single 
criterion can model them all. The task remains, therefore, 
to identify the appropriate failure criterion for each 
failure mode. 

Empirical data are required for the verification of 
any analytical model or failure criterion. A substantial 
data base has been developed characterizing the effects 
of joint geometry [ 9 ] , laminate configuration [9-13] , 
stacking sequence [14], elevated temperature [15,16] and 
moisture concentration [17,18] on composite bolted joint 
performance. Results from these experimental investigations 
provided valuable insight into the failure behavior of 
composite bolted joints which the authors used in the 
formulation of analytical failure models herein. 
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In this study parallr' experimental and analytical 
characterization of bolted joint behavior was conducted 
in an attempt to develop semi-empirical strength models 
for each of the three primary failure modes: net tension, 
shearout, and bearing. A two-dimensional, plane-stress, 
orthotropic, finite-element analysis was used to determine 
the state of stress around the loaded hole in a finite- 
width, semi-infinite orthotropic strip. The stress 
information obtained by the finite element analysis was 
then combined with the failure criteria to predict 
strengths. Observed similarities between the failure 
behavior of notched composites and mechanically fastened 
composite joints suggested that the point stress failure 
criterion might be successfully applied to the shearout 
and not tension failure modes. The bearing failure 
mode has been shown to be dependent on the multiaxial 
state of stress and thus the Tsai-Wu criterion was 
considered. 

It is important to point out that the failure 
criteria proposed in the development of the bolted joint 
strength models are based on intrinsic material properties 
which are experimentally determined using appropriate test 
methods. This leads to the desirable result that the 
design of bolted joints can be accomplished using specific 
material properties instead of empirical bolted joint 
strength data. Thus the development of a reliable, 
comprehensive bolted joir.t strength analysis could lead 
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to substantial reductions in the size of the empirical 
data base required for bolted joint design. 


EXPERIMENTAL PROCEDURE 


* 


Data characterizing the effect of joint geometry 

(i.e. e/D and w/D) on strength and failure mode has been 

generated previously [15] . The experimental program 

conducted in this study was designed to determine the 

T C T 

basic material properties (E x , Ey . v X y , G X y, S x , S x , Sy 
and S„„) of the designated AS/3501-6 graphite/epoxy 
laminate and to determine the notched strength behavior 
for a range of circular notch sizes loaded under 
bolted joint conditions. A set of specimens was designed 
to exhibit each of the three primary failure modes holding 
e/D and w/D constant, while varying bolt diameter, D. 
Several specimens were instrumented for verification of 
the finite element analysis. The details of this 
experimental program are discussed below. 

Specimen Fabrication 

The test specimens for both the material 
characterization study and the fastener size effect 
study were fabricated from a single lot of Hercules 
AS/3501-6 graphite/epoxy prepreg. The laminate 
configuration was (45/0/-45/02/“45/0/45/C>2/9Gl g which 
resulted in a cured panel thickness of 3.0 mm (0.118 in.) 
using manufacturer recommended autoclave processing 
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procedures. All laminates were ultrasonically i'* i pec ted 
for processing related flaws to insure the exclusion of 
the flawed material from the test program. 

The test specimens used for the material 
characterization study were of standard configurations 
for tensile tests, IITRI compression tests and two-rail, 
shear tests. Details of the test procedure, fixturing and 
specimen geometry can be found in the literature [19,20,21]. 
The bolted joint specimen was of the same single fastener 
design described by Wilson and Pipes in an earlier report 
[15] and is shown in Figure 1. 

Test coupons were machined to size using a diamond 
saw to minimize machining-induced damaged. Similarly, 
diamond core drills were used to machine the 3.2 mm 
(0.125 in.), 6.4 mm (0.250 in.) and 9.5 mm (0.375 in.) 
holes in the specimens. The integrity of the machined 
specimens was verified using ultrasonic inspection and 
any flawed coupons were discarded from the test program. 
Width, thickness hole diameter and hole location 
dimensions were recorded for each specimen. 

Testing Procedure 

The general test procedures used in the material 
characterization tests are well documented and will not be 
reiterated in this report. However, specific procedural 
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aspects of the test program are reported herein for 
completeness. 

The specimens were loaded by a floor model TTC 
Instron test machine at a crosshead speed of 0.02 cm/min. 
Strengths were determined by observing the maximum load 
carried by the specimen * :ermined by the load- 

deflection curve. Deformations were measured using 
electrical-resistance foil strain gages and digital strain 
instrumentation. Modulus measurement utilized gages with 
a 6.35 mm (0.1250 in.) grid and point strain measurements 
used gages with a 0.38 mm (0.015 in.) grid length. 

Strain measurements were typically taken at 98 Nt load 
increments. A minimum of three specimens were tested 
to determine each material property and three replicates 
of each bolted joint configuration were tested. 

The testing procedure used for the bolted joint 
tests was similar to that used for the standard static 
tensile property characterization. A specially designed 
clevis fixture (Fig. 2) was used to simulate bolted load 
reaction through the hole in the coupon while standard 
wedge action friction grips introduced load at the tabbed 
end of the specimen. For the pinned tests the load was 
reacted by a high strength steel pin and no constraining 
contact was allowed on the laminate surface surrounding 
the pin. To simulate out of plane constraint annular 
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washers with an inside diameter equal to the pin diameter 
and an outside diameter equal to twice the pin diameter 
were employed. Various sizes of pins and ushers were 
available to insure excellent fit of the washer to the 
pin and the pin to the fastener hole. The torques applied 
to the simulated fasteners were defined as finger tight. 

The effect of fastener fit and torque on bolted joint 
strength is significant and is discussed in later sections. 

Failure Strength Analysis 

The ability of analytical models to describe 
composite bolted joint strength and failure mode 
behavior is dependent upon the knowledge of the state of 
stress and application of appropriate failure criteria. 

A finite-element model was developed to analyze the 
state of stress for the bolted joint geometries and loadings 
considered in this program. Stress profiles predicted by 
the finite-element model along the critical failure 
planes defined in Figure 3 for the bearing, shear out and 
net tension modes were determined and verified by experiment. 
Functional representations of these stress profiles 
including the parametric effects of geometry were then 
developed for combination with the failure criteria. 

Based on simularities observed between notched strength 
and bolted joint failure behavior the effectiveness of the 



"point stress" criterion was investigated for the net tension 
and shearout failure modes. The Tsai-Wu failure criterion 
was studied for the bearing failure mode because of the 
known effects of the multiaxial state of stress in the 
bearing failure zone. Semi -empirically based, closed-form 
relationships result from the introduction of the approximate 
stress functions into the failure criteria which can be 
used to predict composite bolted joint strengths. 

The Structural Analysis Program SAPV was used to 
develop a two-dimensional linear-elastic finite element 
model of a composite bolted joint. The model employed 
the plane-stress, quadrilateral element with orthotropic 
material properties to form the mesh subjected to the 
boundary conditions depicted in Figure 4. The bolted 
joint was treated as a loaded hole containing a frictionless 
rigid inclusion in a finite-width semi-infinite strip. 

Sixteen different geometric variations of the model were 
developed to analyze the effects of geometric parameters 
e/D and w/D on the state of stress. The experimentally 
determined material properties used in the finite element 
model and in the failure criteria are listed in Table 1. 

The stress profiles developed along the three 
critical failure planes are plotted in non-dimensional form 
in Figures 5-7. Only the primary stress component important 



Figure 4. Finite Element Model Used to Analyze the Bolted Joint 
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Table 1. 

Material and Strength Properties 


Property 


Value 



v 

v 

G 

a 


xy 

yx 

xy 

ult 
x , tens 


ult 
y , tens 


a ult 
x , comp 


T 


Ult 

xy 


10.9 10 6 psi 

75.2 

GPa 

4.39 10 6 psi 

30.3 

GPa 

0.397 



0.160 



2.00 10 6 psi 

13.8 

GPa 

131.0 ksi 

903.2 

MPa 

45.2 ksi 

311.6 

MPa 

140.0 ksi 

965.0 

MPa 


21.7 ksi 


149.6 MPa 
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to each failure mode is presented (i.e. shear stress 
for shearout, tensile stress for net tension, etc.). 

To develop confidence in the stresses predicted by the 
finite-element model along the critical planes for the net 
tension, shearout and bearing failure modes, specimens were 
instrumented with strain gages along these planes. The 
gage grid layers were 0.381 mm (0.015 in.) and thus closely 
approximated point strain measurements for comparison with 
the finite-element predicted values. The data plotted in 
Figures 5-7 summarize the results for the net tension, 
shearout and bearing modes respectively. The strain 
results showed relatively poor agreement with the finite- 
element predicted strains along the net tension plane 
(Fig. 5) but a plot of the strains along a plane 0.25 mm 
(0.01 in.) above the net section plane showed excellent 
agreement with the experimental results. Normal strain 
component measurements taken at 45° to the shearout plane 
agree well with the finite-element predictions as shown 
in Fig. 6. The bearing strains predicted by finite 
elements exhibited the proper trend but did not agree closely 
with the measured values (Fig. 7) . The relatively close 
agreement away from the hole tends to indicate that the 
stresses are much greater near the hole than predicted 
by the model. This result could be due to some outer ply 
buckling behavior or inaccurate modeling of the bolt- joint 
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Figure 5. 


Plot of the Strain Profile for a. Along 
the Net Tension Plane ^ 






Based on observed experimental characterization 
of bolted joint behavior a "^oint stress” failure criterion 
[22,23,24] was applied to the net tension and shearout 
failure modes while the Tsai-Wu criterion [25] was 
used for the bearing mode. The "point stress" criterion 
developed by Whitney for tensile strength of a notched 
laminate is based on the stress distribution adjacent 
to the notch. The criterion states that failure occurs 
when the stress at some critical distance, "d 0 ", from 
the edge of the notch reaches the unnotched laminate 
strength. The failure criterion as applied to the net 
tension and shearout failure regimes in bolted joints 
is illustrated graphically in Figures 8 and 9. Only the 
stress component primarily responsible for failure in the 
particular mode under consideration is used in the "point 
stress" criterion. Thus the tensile stress distribution 
0 y(x) across the net tension plane of the fastener is 
applied in the net tension failure mode while the shear 
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of the Laminate 


Figure 8. 


Description of the Net Tension Point 
Stress Failure Criterion 


= Unnotched Shear 
Strength of Laminate 


Figure 9. Description of the Shearout Point Stress 
Failure Criterion 
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stress distribution, x (y) along the shearout plane 

xy 

is used for the shearout mode. 


The Tsai-Wu failure criterion may be expressed 
in the following form to predict strength of composites 
subjected to multiaxial states of stress, 

p i°i + ± 1 

where and F^j are strength tensors of 2nd and 4th rank, 
and the are usual contracted notation for stress 
components. The strength tensors must be determined by 
experiment for eact. material system considered. 
Specifically, the strength tensors are defined as 
follows . 

F, = — + — 

1 x t x c 
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x t x c 


F = i + JL 

2 Y t Y c 


22 


Y t Y c 


F 6 = 0 


66 


1 

7 


where 

x fc = laminate tensile strength in the load direct. 

x c = laminate compressive strength in the load direct. 

Y. = laminate tensile strength transverse to the 
loading axis. 
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Y * laminate compressive strength transverse to 
^ the loading axis 

S * inplane shear strength of the laminate 

The quantity cannot be measured by uniaxial tensile 
tests but can be measured by a biaxial test. Failure 
occurs when the sum of the tensor products on the left 
side of equation 1 equals or exceeds unity. 

Implicit in the application of either of the 
above failure criteria is an analytic representation of 
the stresses in the joint. Since a closed form solution 
to the problem is not known, finite element analysis was 
used to determine the stress profiles. 

A simple collocation method was used to fit the 
stress distributions by polynomial expansions of the space 
variable in nondimensional form. Net tension stress 
distributions were fit by a polynomial of the following 
form: 


— = A + B 

0 


f?.l 2 + c f— I 

(yj [yJ 


+ D 


+ E 


8 


( 2 ) 


where the constants A through E were determined by 
satisfying equilibrium and from points on the finite 
element stress distribution curve (Fig. 10) . Agreement 
between the finite-element results and the approximate 
functional representation generated by equation (2) is 
shown in Figure 10. The shearout stress distribution 
was modeled by an expression of the form: 
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Figure 10. Finite Element Predicted Net Tension Stress 
Distribution and Comparison to Approximate 
Polynomial Fit 
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+ c 

o 

X 

k 4 


where the constants A, B and C were found by fitting 
three data points at critical locations on the shear 
stress distribution curve, Fig. 11. Comparision of 
the finite element curve with the polynomial 
approximation shows good agreement in the critical region 
of high stress concentration near the hole with much 
poorer agreement in the far field region. The inclusion 
of higher order terms in the polynomial reunited in the 
generation of harmonics which induced significant 
fluctuations in the stresses for small x (near the hole) . 
Only slight increases in accuracy in the functional 
representation of the stress profile were achieved by the 
use of more complex functions, and therefore the 3rd 
order polynomial was considered adequate. 

Parametric studies of the net tension and shear 
stress distributions as functions of the geometric 
parameters e/0 and w/D resulted in the sets of curves 
in Figures 12 through 14. There was no effect by e/D 
on the net tension stress distribution while 
Fig. 12 shows a significant effect of w/D. The 
shape of the curve describing the stress distribution 
did not change but was shifted downward with decreasing 
w/D. Similar effects were noticed for the shear stress 
distributions. Increasing e/D resulted in a shifting of 



to Approximate Polynomial 










Figure 14. The Effect of Varying W/D on Shearout Stress Profiles for e/D 
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the curve downward without changing its shape or shifting 
the location of the maximum stress (Fig. 13) . Increases 
in w/D caused a linear increase in the magnitude of the 
stresses without changing the shape or location of the 
curve (Fig. 14). These results for the net tension and 
shearout stress distributions are better seen by plotting 
the peak stress versus e/D and w/D as shown in Figures 
15 through 17. Non-dimensional parameters were developed 
from the relationships describing these geometric effects 
on the stress profiles and used to incorporate the 
geometric parameters w/D and e/D into the functional 
stress profile representations. Defining a function 
K (w/D) and fitting the curve in Fig. 15 it is found that 

Uw/D) » 0.956 - 0.044 x + 0.089 x 2 
where x = w/2D 

for the net tension mode. Similarly it can be shown that 
two parameters p (w/D) and C (e/D) can be defined and found 
to have the following forms derived from the curves in 
Figures 16 and 17 

p ,„ /0 ) - 

T max 

w/D**6 

where M = slope of the t „„ versus w/D curve in 

max 

Figure 16 and 

C (e/D) = 0.75 e (0,5 " e/D) + 0.84 
for the shearout failure mode. 
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Figure 17. The Variation of Maximum Shear Stress Concentration with e/D 
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Combining the geometric parameters and the 
original stress profiles results in functions which 
describe the stress profile incorporating the 
geometric effects. 



i (w/D) 


A + 


B(|) 


C(|) 


D(|) 


E(|) 


and 


sa 


* C(e/D) *p (w/D) 


x n x^ 

a( tt» + B <-r> + c( ir> 


These approximate expressions for the stress are suitable 
for incorporation into the failure criterion. 


The "point stress" 
expressed as 


x“d„ 


failure criterion is 


where 


o is the stress component adjacent to the hole 
for the failure mode under consideration 

o ) is the unmatched strength of the laminate for 
the failure mode under consideration. 

d^ is the critical distance, 
o 

The critical distance parameter is a function 
of hole size and has the following form (22] . 



Constants C and m are found from empirical data by finding d Q 
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for two notch sizes. This was carried out for both failure 
modes with different resultant m and C values. The values 
of m and C for the net tension failure mode were 0.83 
and 8.30, respectively, while for the shearout failure 
mode m was 0.140 and C was 7.12. 

The bearing mode failure criterion considers the 

effects of all in-plane stress components in the compressive 

region of the joint. Each stress component a , a and 

x y 

T x y are plotted in Figure 18 from the centerline 
of the hole's edge along the y-axis. It is apparent 
that the maximum compressive stress o x does not occur 
at the centerline but very slightly off the centerline. 

This must be considered when applying the Tsai-Wu failure 
criterion. 

Parametric studies evaluating the effects of 
geometry on the bearing stress profiles are summarized 
in Figures 19-22. The plane stress finite element 
analysis shows that increasing w/D results in increasing 
compressive stress and a change in the distribution of 
the stress (Fig. 19) . Similarly Fig. 20 shows that 
increasing e/D results in increasing stresses and slight 
changes in their distribution. This same information 
is represented in a more convenient form by Figures 21 
and 22. Notice that the maximum predicted compressive 
stress is less than half that of the stress calculated 
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for the projected pin bearing area. This leads to further 
suspicion that two-dimensional, plane-stress analysis 
does not accurately predict the stresses in the bearing 
region. 


cr/cr. 
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2.0 
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' 0.0 0.4 0.8 1.2 1.6 2.0 

Distance from Pin x/D 

Figure 20. Compressive Stress Distribution Along 
Joint Centerline as a Function of e/D 







RESULTS AND DISCUSSION 


The experimental data from the fastener, hole 
size study are reported in tabular form in Appendix A. 

For the limited range of fastener sizes investigated 
the effect was small. The net tension failure exhibited 
slightly decreasing strength with increasing hole diameter. 
Shearout failures resulted in the opposite trend, 
increasing strength with increasing hole size. Bearing 
failures showed no hole size effect at all. 

The hole size effects for the net tension and 
shearout modes were incorporated into the failure model 
through d Q . In Figures 23 and 24 failure curves were 
calculated using the derived failure criteria for a 
range of fastener sizes while holding w/D and e/D constant. 
The excellent correlation found between the experimental 
data and the analytical results for both failure modes 
occurs because two of the three data points were used in 
the determination of constants in the failure model and 
thus is not surprising. Since the hole size effect is 
small, the predicted behavior should remain accurate 
for the range of fastener sizes of practical interest. 
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Figure 23. Correlation of Model Predicted Net 
Tension Strengths with Experimental 
Data as a Function of Hole Diameter 






r 0 - Unnotched Sheor Strength 

T so“ Par Field Stress 

«■■■* Theoretical Model 
m*O.II4 c*-U3 
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Prom previous test programs [15,9] data was 
available to investigate the accuracy of the newly developed 
failure criteria in predicting net tension and shearout 
failure strength for several w/D geometries. In 
Figure 25 the solid curves are the predicted failure 
surfaces for two distinct fastener sizes for the net 
tension mode. The circular symbols represent data from 
the current test program. The triangular symbols are 
data from Hartsmith [2] for a similar laminate but from 
a different graphite-epoxy system (Narmo 5208/T300) . 

The agreement between the experimental data and the 
predicted failures is excellent. Similar results were 
found for the shearout failure mode as seen in Figure 26. 

Bearing failure analysis was attempted using the 
Tsai-Wu strength criterion. For convenience the joint 
geometry was redefined in polar coordinates as shown in 
Figure 27. The Tsai-Wu criterion was applied to the 
laminate at each element location on the hole boundary 
between the angular positions 0*0° and 90°. The 
state of stress at failure was determined by setting the 
finite-element far field stress input equal to the far 
field stress experiment for bearing failure. This set 
of stress components obtained were then used as input 
to the Tsai-Wu criterion to test its validity for bearing 
failure analysis. 



Figure 25. Correlation of Net Tension Strength as a Function of W/D with 
Experimental Data 



f , 1 


Figure 26. Correlation of Shearout Strength as a 
Function of W/D with Experimental Data 
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The Tsai-Wu criterion predicted strength was seen 
to increase monotonically with increasing angular position 
in the 0° to 45° range while developing an exponential 
character and increasing without bound as the 90° 
position was approached. This behavior presented a 
problem in interpretation of results since the bearing 
load reacted by the laminate decreases as the 90° position 
is approached. Since the 0° to 45° interval reacts most 
of the bearing load which causes failure ahead of the 
fastener this region was selected as the valid region for 
application of the criterion. The maximum value of the 
criterion in this region was 0.5 at the 0°±45° location 
when the "failure stress components" were input. This 
value is significantly different from the value of 
approximately 1.0 expected at failure and is non-conservative. 

One possible explanation for the poor results can 
be linked to the plane stress finite element analysis. 

The plane stress assumption for the bearing analysis is 
not a good approximation. This was verified for the pin 
bearing case where measured strains were much greater 
than predicted strains. The interlaminar stress components 
°z ' T xz ' an< * T yz cou - !Lc * significantly affect the compressive 
strength of the laminate. Another factor which affects the 
results is fastener torque and washer size. According to 
Stockdale and Matthews [26] in Figure 28 a twofold increase 
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(Stockdale and Matthews [26]. 
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in strength can be achieved by increasing the fastener 
clamping load from finger light to approximately 14 MN. 
Similarly if the washer diameter is increased from 2 to 
4 times the fastener diameter (D) the failure strength 
exhibits a twofold increase. Thus it is possible that 
a combination of low torque and small washer area 
produced low joint bearing strength. A check with the 
literature appears to refute this hypothesis, however. 

The bearing strengths which were measured in the present 
tests compare closely with those of Hartsmith and others 
[9,10]. With the accuracy of the bearing stress analysis 
questionable it is not possible to assess the validity of 
the Tsai-Wu criterion for bearing analysis. 

The models developed for the net tension and 
shearout failure did not deal with changes in material 
anisotropy and the resultant effects on the strength. 

If the anisotropy (orthotropy) is defined by the ratio 
of inplane longitudinal modulus to transverse modulus, 
the effects on the *nplane stress distributions can be 
studied as a function of anisotropy ratio. Some preliminary 
work indicates that the primary effect of anisotropy is to 
change the shape o e the stress distributions. The location 
of the point of maximum stress remains the same however. 

By formulating approximate stress functions for the 
different anisotropy conditions the models can be 
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incorporated into the failure criteria without changing 
the accuracy of the strength predictions. More research 
both analytical and experimental needs to be done in 
order to fully understand the effects of anisotropy. 


CONCLUSIONS 


The failure behavior of composite bolted joints 
has been investigated for the three primary failure 
modes, bearing, shearout and net tension. Strength 
analyses were developed for each mode and their accuracy 
determined by comparison with empirical data. 

The analysis of composite bolted joint strength 
is dependent upon determination of the state of stress 
in the joint and the application of an appropriate 
failure criterion. A two-dimensional, plane stress, 
finite-element model was found to accurately predict 
stresses along the net tension and shearout planes as 
verified by comparison with experimental results. The 
plane stress assumption appears to be invalid for analyses 
of stresses near fastener hole in the bearing region of 
the joint. The predicted stresses appear to be approximately 
one half the magnitude measured by strain gages. Since 
the interlaminar stresses could significantly affect 
the strength of the laminate in this region it is suggested 
that a three-dimensional analysis be developed. 

Parametric studies performed to evaluate the 
effect of joint geometry on the state of stress around 
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the joint revealed the existence of relatively simple 
relationships between maximum stress and geometry. This 
allowed development of functional form approximating the 
stress profiles fcr each failure mode. This development 
of analytical approximation of the stresses as a 
function of joint is important in the reduction of the 
empirical data base needed in composite joint design. 

Semi-empirically based strength analyses were 
developed for each failure mode by application of the 
failure criteria and the stress analysis. After 
comparison with experimental data it was concluded 
that failure models using the "point stress" criterion 
for net tension and shearout failure modes accurately predict 
bolted joint strengths for the laminate configuration and 
geometries investigated. The questionable accuracy of 
the plane stress analysis for the bearing made prohibited 
any conclusions about the applicability of the Tsai-Wu 
failure criterion. 

The development of an approximate closed form 
strength analysis based on the "point stress" criterion and 
incorporated geometric effects results in a potential 
reduction in the empirical data base required for bolted 
joint design. The expansion of this concept to include the 
bearing mode and effects of material anisotropy and 
environment could lead to the development of a comprehensive 
bolted joint design procedure. 
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